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We analyze the appearance of non-Markovian effects in the dynamics of a bipartite system cou- 
pled to a reservoir, which can be described within a class of non-Markovian equations given by a 
generalized Lindblad structure. A novel master equation, which we term quantum Bloch-Boltzmann 
equation, is derived, describing both motional and internal states of a test particle in a quantum 
framework. When due to the preparation of the system or to decoherence effects one of the two 
degrees of freedom is amenable to a classical treatment and not resolved in the final measurement, 
though relevant for the interaction with the reservoir, non-Markovian behaviors such as stretched 
\ exponential or power law decay of coherences can be put into evidence. 

^ ■ PACS numbers: 03.65.Yz, 05.20.Dd, 03.75.-b, 42.50.Lc 
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The complete isolation of quantum mechanical systems, which should arise by a perfect shielding from the envi- 
ronment, can of course in general only be an idealization. The study of the dynamics of open quantum systems then 
\ naturally becomes of great interest [H, especially when it comes to a realistic description of experimental situations. 
II While for the case of a closed quantum system the time evolution is given by a one-parameter unitary group charac- 
terized by a self-adjoint Hamiltonian, the situation is more involved for an open quantum system, where dynamical 
. evolutions including irreversible effects like dissipation and decoherence must also be considered. The possible struc- 
^ ■ tures of dynamical equations for an open quantum system are not known in full generality, despite the huge efforts 
O^, devoted to the problem. A well-known result of paramount importance has been obtained for Markovian dynamics, 
asking the mapping giving the dynamics to be a completely positive quantum dynamical semigroup. The expression 
\ of the generators of such semigroups, which gives the master equation for the statistical operator of the system, has 
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I. INTRODUCTION 



been fully characterized 0, 0], providing a reference structure, often referred to as Lindblad equation. Such Lindblad 
type master equations ensure a well-defined time evolution, preserving in particular the positivity of the statistical 
\^ , operator. The different terms and operators appearing in it are often naturally amenable to a direct physical interpre- 
' tation. Moreover analytical approaches are often feasible, and when this is not the case numerical studies can always 
, be performed, by considering Monte Carlo simulations of suitable stochastic differential equations associated to the 

master equation via a particular unraveling. 
\ Such a general and physically transparent characterization is not available for master equations describing a non- 
25 , Markovian dynamics. However systems exhibiting non-Markovian dynamics, such as memory effects and decay behav- 
iors other than exponential, are also of great interest both for practical applications and from a conceptual standpoint. 
In this spirit major efforts have been devoted to derive possibly general classes of master equations, which while pro- 
viding well-defined time evolutions also describe non-Markovian effects. Various difhculties appear in this connection. 
^ j , In particular it is important to provide a link between the operators entering the structure of generalized master 
' equations and quantities of physical relevance characterizing the environment and its coupling to the system. General 
classes of non-Markovian master equations have been obtained in the literature [1, [E, IB] j also pointing to possible 
physical applications. In particular the analysis of the interaction of a quantum system with a structured reservoir, 
performed via a time-convolutionless projection operator technique relying on the use of correlated projection opera- 
tors adapted to the structured reservoir Q , has led to point out a generalized Lindblad structure Q . This generalized 
Lindblad structure describes a non-Markovian dynamics on states which are given by classical convex mixtures of 
subcoUections, that is positive trace class operators with trace equal or less than one, naturally appearing in the 
description of quantum experiments @. Master equations of this form have already been proposed in an utt erly 
different context in order to introduce the notion of event in the description of quantum mechanical systems [lO|, 
for the purpose of better understanding the interplay between classical and quantum description of physical reality. 
More recently and to the point similar equations have been considered for the statistical operator of an active atom 



'Electronic address: ' bassano .vacchiniOmi . infn . it I 



2 



interacting through colhsions with a gas, when describing in a classical way the centre of mass degrees of freedom 
. . . 

General physical mechanisms leading to the appearance of such generalized Lindblad structures which can account 
for non-Markovian effects have already been conceived. This is the case if one studies the dynamics of an open 
system coupled to a structured reservoir using the abovementioned time-convolutionless projection operator technique, 
provided the projectors used in obtaining the reduced equations of motion do project on classically correlated states 
between system and environment, rather than simply on a factorized state, as in the common wisdom 0, [T^ [l3| . 
Another natural situation leading to this generalized Lindblad structure appears in what has been called generalized 
Born-Markov approximation Here one considers the usual second order perturbation scheme, but once again the 
state of system and bath is supposed not to be factorized, but rather given by a convex mixture of factorized states. 
The indexes of the mixture are related to the structure of the bath and of the interaction Hamiltonian. Further work 
has traced back the derivation of non-Markovian equations of this form to the existence of extra unobserved degrees 
of freedom mediating the entanglement between the considered system and a Markovian reservoir [l6| . Earlier 
work also led to this kind of non-Markovian master equations for a system interacting with an environment with 
a finite heat capacity, so that energy exchanges between system and reservoir also affect the energy distribution of 
the reservoir. It has been recently shown that the same quantum master equation can also be derived in a physically 
more transparent manner by means of the projection superoperator technique (l8| . 

In the present paper we show how such generalized master equations naturally arise by considering a bipartite 
quantum system interacting with a Markovian reservoir, whenever decoherence effects or superselection rules affect 
only one kind of degrees of freedom of the bipartite system. A non-Markovian behavior then appears when the 
thus emerged classical label is not resolved in the final measurement. The analysis is done by means of a concrete 
and relevant physical example. We consider the dynamics of a quantum test particle, whose internal and centre of 
mass degrees of freedom are both described quantum mechanically, interacting e.g. with a gaseous background. The 
appearance of non-Markovian features is related to the involvement of both internal and centre of mass degrees of 
freedom in the scattering amplitude which describes the coupling between bipartite system and environment. As a 
first step a novel quantum master equation is heuristically derived, which extends previous work on the quantum linear 
Boltzmann equation [1^, |23, [2l|, [2 3 . focussin g on a quantum description of the centre of mass degrees of freedom, and 
on the Bloch-Boltzmann equation [ill [23[ . [2^ l25l , which describes in a classical way the motion of the test particle, 
but retains a quantum expression for the dynamics of its internal degrees of freedom. According to this terminology 
the obtained equation is termed quantum Bloch-Boltzmann equation. Two limiting situations then naturally appears. 
When decoherence affects more strongly or equivalently on a shorter time scale the centre of mass degrees of freedom, 
a generalized Lindblad structure corresponding to the Bloch-Boltzmann equation appears, in which the momentum of 
the test particle is treated classically. This equation describes non-Markovian effects when the final measurement only 
affects the internal degrees of freedom. In a similar way, when the experimental effort is devoted to study quantum 
superpositions of motional degrees of freedom, e.g. in interferometers for massive particles studying robustness of 
their quantum behavior, non-Markovian features can appear if the internal degrees of freedom influence the collisional 
scattering cross section but are later not observed in the assessment of the interference pattern. 

The paper is organized as follows. In Sect. [TT|we describe mathematical framework and expression of the generalized 
Lindblad structure, in Sect. IIIII we outline the derivation of the quantum Bloch-Boltzmann equation, involving both 
centre of mass and internal degrees of freedom, starting from the classical linear Boltzmann equation and the known 
expression of the quantum linear Boltzmann equation involving only the motional degrees of freedom. Sect. IIVI is 
then devoted to consider the reduced dynamics of either internal or centre of mass degrees of freedom, for which a 
generalized Lindblad structure follows, showing by means of example the appearance of non-Markovian behaviors. 
Finally in Sect. [V|we briefly comment on our results. 

II. GENERALIZED LINDBLAD STRUCTURE 

We now want to introduce the abovementioned non-Markovian generalization of the Lindblad structure, which 
is easily obtained considering the standard theorem of Gorini, Kossakowski, Sudarshan and Lindblad for a special 
choice of Hilbert space for the open system and of the expression of its statistical operator. Let us consider a bipartite 
quantum system described on a Hilbert space 7i (Xi Tis, with Ti. and Tis separable. Exploiting the isomorphism of 
Hb with either of the possible physically relevant choices of Hilbert space, such as (D", P ((D) or (R), the tensor 
product can be expressed as a direct sum or direct integral, thus naturally introducing a label a. For the case in 
which Hb describes a system with n degrees of freedom one can consider the two equivalent constructions of the same 
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bipartite Hilbert space 

n 

H®C" = 07^, (1) 

a = l 

and similarly for P (C), replacing the finite sum with a series. On similar grounds a continuous index appears for the 



case 



T-L' 



(M) = / dan, (2) 



exploiting the notion of direct integral of Hilbert spaces (see e.g. [23]). It is then of course possible to consider a 
statistical operator for the bipartite system, that is to say a positive, trace class operator on 7i (8" Tis, normalized to 
one. Considering for the sake of example the case oiTi ® (D" = ^ind denoting by g the statistical operator 

of the system one can consider the general expression of the Lindblad equation for g. Let us now restrict however 
to statistical operators whose matrix representation is block diagonal, so that they can equivalently be written as 
g — X]a=i Pa®\c»){a\ or as p = {pi, . . . , pa, ■ ■ ■ , Pn)- The index a can now really be interpreted as a classical label 
indexing the various subcoUections pa € TC(7i), which are given by positive, trace class operators on 7i with trace 
less or equal than one. Such a block diagonal statistical operator 

g = (Pi, . . . . . . 

fixed by the set of subcoUections {pa}^^! n normalized according to 

n 

Trp = ^Tr-HPa = l- (3) 

Q = l 

The set of trace class operators which are block diagonal is a closed subalgebra of the set of all trace class operators, 
whose dual space is given by the closed subalgebra of bounded operators also having a block diagonal structure. 
Equivalently one can say that this subclass of statistical operators only provides information on the expectation 
values of observables diagonal with respect to the label a, which thus correspond to the only relevant variables for 
the system under consideration. Considering a statistical operator in the subalgebra of block diagonal trace class 
operators p G TCdiag ® (D") and an observable given by a block diagonal bounded operator B g ;Sdiag C") 

B = (Bi, . . . , Bq, . . . , B„) , 

with Ba E B (7Y), the duality relation is given by 

n 

{B,g) = ^Tr„(B„p„). (4) 

It is now of interest to consider the expression of the generator of a quantum dynamical semigroup acting on 
this bipartite space when applied to such block diagonal states or observables, with the further constraint that the 
time evolved state or observable still preserves this simple block diagonal structure, thus defining a dynamics which 
remains within the spaces TCdiag CH (8) C") or Sdiag (8) C") respectively. In the Schrodinger picture this generalized 
Lindblad structure can be written in terms of coupled equations for the different subcoUections Pa according to ^ 



dt^" " h 



d i 

A 13=1 



(5) 



leading due to the duality relation Eq. ^ to the following equations in Heisenberg picture for the components Bq of 
a block diagonal observable 



Ab„ = +i[H",BJ+EE[RrB,Rf -i{RftRf,B4 



A /3=1 



(6) 



In Eqs. ^ and ^ the index a runs from 1 to n, the operators H" are self-adjoint on H and R^" are operators on 
7i, with A a further index labelling the various Lindblad operators. 
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Introducing the mapping 



. , d d d 

= (d^^^"-"d^^----'dt^" 



one can therefore write for the time evolution in Schrodinger picture 

= e*^e(0) 

= e*^(pi(0),...,pa(0),...,p„(0)), 

and similarly for the Heisenberg picture using the mapping C dual to C according to the relation Eq. ([4]). Now Eq. 
([5]) provides a Markovian set of equations for the statistical operator Q{t) = {pi{t), . . . , pa{t), . . . , Pn{t)) on H (Xi (D", 
but a non-Markovian dynamics for the statistical operator 

n 

which is a statistical operator on the Hilbert space Ti. only. In particular it is not possible to define a mapping from 
w{0) to wit) according to the non-commutativity of the following diagram 

giO) = (pi(0), . . . , p„(0))'-^ g{t) - . . . , Pnit)) 



^(0) = ELi Pa(o) ^ ^(t) - ELi W 

arising because of the loss of information in going from g(t) to w{t). The set of equations given by Eq. ([S]) thus 
provides a non-Markovian dynamics for the statistical operator wit) supposed to be expressible at any time as a 
mixture of subcoUections pa (t) , or equivalently as a convex combination with weights Pa (t) = Tt-h pa (t) of statistical 
operators given by Wa {t) — Pa (t) / Tr-H Pa {t) . This last standpoint stresses the appearance of the classical probab ility 
distribution {Pa{t)}a=i „i which justifies the name random Lindblad equations or Lindblad rate equations (l^.ll6j|. 
also given to equations falling within the class given by Eq. ([5]). 

The statistical operator w{t) can arise in a twofold way. Either by taking the trace of a block diagonal g{t) with 
respect to (D", corresponding to a situation in which one considers a reduced dynamics of the bipartite system with 
respect to the degrees of freedom which behave effectively in a classical way, or by assuming that the state of the 
system under study living in the Hilbert space Ti. is specified at the initial time as a convex combination of n statistical 
operators with suitable weights, and retains this form throughout the dynamics. The first type of realization makes 
it intuitively clear why Eq. ([5]) encompasses non-Markovian situations. By looking at the time evolution of w(t) 
only one is considering a restricted set of variables with respect to the full collection {pa {t)}a=i,...,7i: for which the 
time evolution law would be Markovian. The set of relevant physical variables then determines whether or not the 
dynamics is Markovian. Statistical operators of the form given by Eq. ([7]) naturally appear in connection with a 
structured reservoir, the label a being then connected to a characterization of the reservoir itself, e.g. labelling different 
energy bands. This is the case both when considering a projection operator technique assuming classical correlated 
states between system and reservoir [7] , or more simply in the so called generalized Born-Markov approximation [3] j 
using a classically correlated state in the derivation of the master equation to second order in the perturbation. As 
we shall argue below statistical operators of this form also appear when considering a bipartite system interacting 
with a reservoir, when due to decoherence one of the two kind of degrees of freedom behaves classically, and despite 
characterizing the initial preparation and being relevant for the interaction with the environment, cannot later be 
resolved by the measurement apparatus. Needless to say the formal scheme developed above can also be implemented 
when the label a runs over a countable set or even in the continuum, corresponding to a Hilbert space construction 
as depicted in Eq. 



III. DERIVATION OF QUANTUM BLOCH-BOLTZMANN EQUATION 



We now address the issue of the derivation of a master equation describing the dynamics of both internal and 
centre of mass degrees of freedom of a quantum test particle interacting through collisions with a reservoir such as a 
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background gas. Here we will outline a heuristic derivation, leaving a more detailed and microscopic one for a later 
publication, in particular we will concentrate only on the incoherent terms corresponding to gain and loss term in the 
classical case, leaving aside the Hamiltonian contributions due to kinetic term and forward scattering. To do this we 
will build on the known results for the quantum linear Boltzmann equation (lol . [20l [2ll . [2^ , describin g in a quantum 
framework the centre of mass degrees of freedom, and on the so called Bloch-Boltzmann equation pl].l23l. [2^ . [25l . , 
which accounts for a semiclassical description of both internal and centre of mass degrees of freedom. As a starting 
point we take the equations of motion considered in [TTI| for the collection {p{P)}p^^3 of trace class operators on 
the space (D" of internal degrees of freedom, labelled with the continuous index P now characterizing the classical 
momentum of the test particle. This result extends the work i n [23l[26| by considering also the dynamics of the centre 
of mass degrees of freedom of the test particle, which in [H, |26|| is supposed to have infinite mass and therefore to 
be at rest. It is most convenient to write the equation putting into evidence the constraints due to momentum and 
energy conservation in the single interaction events, thus obtaining 



ijkl 

£ij—£kl 



/ pf2 n p2 2 \ 

\2M 2m ^ 2M 2m J ' 

(rel (p, P) , rel (p', P')) fki (^1 (p, P) , rcl (p', P')) E,,p{P')eIi 



X 

1 ^gas 

2 ml 



JdP' J dp' J dpppip) 



ijkl 



/ p2 2 p/2 12 \ 

X<5 + + fiiJ, ^ huJi]S^(P+p-P' -p') 

\2M 2m ' 2M 2m ) ^ ^ ^' 

xU, (rel (p', P') , rel (p, P)) fn (rel (p', P') , rel (p, P)) [eIe,, , p{P)] . (8) 

In this equation M and m denote the mass of test particle and gas particles respectively, is the reduced mass, 
Hp (p) is the Maxwell-Boltzmann distribution of the gas 

/^/3(P) = 3/2 3 I I ' (9) 

7r3/2p3 \ Pfi) 

with pp — yj2m/ (3 the most probable momentum at temperature T = 1/ (fcB/3), relative momenta are denoted as 



rel (p,P) ^ ^p-^P, (10) 



and 



/y(p,p') = f{p',3^P,i) (11) 

indicates the complex scattering amplitude for a transition from an in state with labels p',j to an out state with 
labels p,i. Finally the matrices Eij correspond to the mappings between energy eigenstates, providing a basis of 
operators in (D", according to 

E^, = (12) 

while hiOj is the energy of the j-lcvcl and £ij = fkui — fkuj denote the possible transition energies. Using the delta of 
momentum conservation Eq. ([H]) can also be expressed using as a variable the momentum transfer Q = P' P in 
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the single collisions, thus coming to 



dt 



ijkl \ / 



ijkl 
£ij—£ki 

X (rel (p, P) , rel (p, P) - Q) f*^ (rel (p, P) , rel (p, P) - Q) E,,p{P + Q)eI 

ijkl \ / 

^ij — £ kl 

X/,, (rel (p, P) - Q, rel (p, P)) (rel (p, P) - Q, rel (p, P)) p(P)} . 

Introducing now by the notation HQ and ± Q the component of a vector parallel and perpendicular to the momentum 
transfer Q, so that P\\q = {P ■ Q) Q/Q^ and P^_Q = P P\\q respectively, one has using Eq. (fTU|) in the delta of 
energy conservation 

^ I PTTT + — — - — + '^y ^ '"'=1 (P||Q' -PlIQ + ' (13) 



2M 2m 2Af 2m / V 2m* m* 

so that in the integral one can use the replacement 



and therefore 



x/„ (rel (p,„.P,«) + f + ^0,„1 (p,«,P,«) - I + 

yJu (.01 (P,«, Px«) + f + (Px,. Px«) - f + 

yEt,f(P + Q)El, 

4^ i: /do/*''.M<|^-£ »'(''.o-M+^" 

X/. (rel (p,^, P.c,) - I + rel (p,^, P,^) + f + 

X At. (rel (p,„ P.^) - f + ^Q, rel (p,^. P.,) + f + 



{i?Li?.„p(P)} 
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One can now perform the translation p ^ p + mQ / (2m*) + m,P\\Q/M + mSjiQ/Q"^ in the gain term, and similarly 
for the loss one, which does not affect the argument of the scattering amplitudes, thus obtaining 



ijkl 



xEijp{P-Q)Eli 

1 rigasm 

2 m2 



X (rel (p,^, - f + Q^Q, rel (p,^. P..,) + f + 

X (rel P.g) - I H- ^Q, rel (p,^, P.g) + f + 

x{eIe,„p{P)}, 

where we also performed the change of variables Q —Q in the gain term and used the simple relation m/ (2m*) — 
m/M = 1 - m/ (2m,). Noting that 

J dpg{p)6{Q-p) = ^J^ dpg{p^Q), (15) 

where the integration on the r.h.s. is restricted to momenta of the gas particle perpendicular to the momentum 
transfer we obtain the equation 



d_ 
dt 



£ij=£kl 

(rel P.g) - f + rel (p,^, P.g) + f + 

xfu (-1 (P.g, P.g) - f + ^Q, rel (p,g, P.g) + f + 



xEijp{P-Q)E, 



kl 



1 ^gas^^ 



7 —^kl 



X /. ( rel (p,^. P.,) - f + ^Q, rel (p,^, P.^) + f + 



X m (rel (P.g, P.g) - f + rel (p,^, P.g) + f + 

x{4,Ei,-,p(P)}. (16) 

Let us now recall the expression of the dynamic structure factor for a gas of free particles obeying Maxwell-Boltzmann 
statistics, which is given by 

[P^l ( (3 {Q^ + 2mEf \ 
Smb{Q,E) = V^Q-P(-^ Q2 J' (17) 



8 



where the variables Q and E denote energy transfer and momentum transfer in a scattering event. The dynamic 
structure factor is a two-point correlation function appearing in the expression of the scattering cross section of a 
probe scattering off a macroscopic sample, in the present case the gas, expressed in terms of momentum and energy 
transferred in the collision '28', ^9*1 . Its general expression is given by the Fourier transform with respect to energy 
and momentum transfer of the density-density correlation function of the medium, and for the case of a sample of non 
interacting particles can be analytically evaluated to give Eq. ()17|) . The physical meaning of the dynamic structure 
factor for the characterization of scattering of a test particle off a gas explains its natural appearance in the expression 
of the quantum linear Boltzmann equation, as already recognized in ,^20, 22, 30, 31]. As we are now going to show 
the dynamic structure factor also appears when considering internal degrees of freedom, the energy transfer being 
now also related to the energy absorbed or released as a consequence of internal transitions. Exploiting Eq. (|15|) we 
observe in fact the identity 

mQ ^ , £^j ^\ _ ( fQ^ + 2m{E{Q,P)+£,j)\Q 



leading via Eq. ^ and Eq. ^ to 

(p^Q+£^f + ^jP\\Q + - M/. (p^q) Smb {Q,E{Q,P)+S.,) , (18) 

where up {p±q) denotes the Maxwell-Boltzmann distribution in two dimensions. The quantity 

E (Q, P) - 2A/ " 2M - 2M + ^iT ' ^^^^ 

actually only depending on P\\q, is the energy transferred to the centre of mass in a collision in which the momentum 
of the test particle changes from P to P + Q. Relying on Eq. (fT5|) and Eq. we can therefore finally write 

-^p(P) = ^ V / dQ / dp^i0{p^Q)SMB{Q,E{Q,P-Q)+8,,) 

Q £kl 

2" QVm 

V Uq [ dpfifj{p^Q)SMBiQ,E{Q,P)+£,,) 

,,,,, J Jq^ 



1^ 


F 




/m. 


Q 


1 Ski 


2 


/m. 



^fki rel (P±Q, P±q) - ^ + T^JT— rcl {p^q,P±q) + - + -t^t^Q 



xE,,p{P-Q)eIi 

1 ^gas 



2 ml 

£ij—£ki 



X/., (rel(p,^,P,p) - f + q^Q'-I(p.Q^^'.q) + | + 



1^ 


F 




/m. 


Q 


1 Ski 


2 


/m. 



-Q,rel(p^Q,P^Q) + ^ + -gi-g 



2 QVm* 

[eIe,,,p{P)}, (20) 



where because of the delta of energy conservation Eq. ()13|) exploited in coming to this final expression the transition 
energy Fij must be equal to zero whenever Q is equal to zero, due to the fact that we consider the gas particles as 
structureless. We stress the dependence on P — Q in the gain term with respect to P in the loss term. Note also 
the very natural appearance of the argument E (Q, P) + Fij in the dynamic structure factor, corresponding to the 
energy transfer in the interaction events, due to both the momentum exchange and the internal transition, whenever 
the scattering is not elastic. 

Eq. (j20p is equivalent to Eq. ([H]), but it is written in a more convenient way for the sake of considering a 
quantum description of the centre of mass degrees of freedom. The quantum master equation for the dynamics of 
both internal and centre of mass degrees of freedom has to be of Lindblad form and to coincide with the semiclassical 
expression (1201) when considering the diagonal matrix elements in the momentum representation. Moreover due to the 
homogeneity of the gas the equation has to reflect the physical invariance under translations, which is expressed at 
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the level of the master equation by the property of covariance under translations, corresponding to the fact that the 
generator of the master equation commutes with the generator of translations. This property has been considered at 
a formal level in [H, [H, [sj, [s^ , leading to a general mathematical characterization of Lindblad structures complying 
with translational invariance, and discussed in a physical framework in (20l . [sij . In view of these requirements the 
quantum master equation is simply obtained by making operator-valued the relevant physical expressions appearing 
in the equation and depending on the momentum of the test particle, such as dynamic structure factor and scattering 
amplitude. In this way one obtains the following master equation for a statistical operator g on the space (R"^) (g)(D", 
which is manifestly in Lindblad form 



dt' 



n 
m 



gas 
2 



ijkl 



/ \ Q 

-QM (p_lq,P-lq) + ^ 



Q"^ /m 



-Q 



Smb (Q, E (Q, P) 

Q 



dQ 



E,,)E,, qeIi ^Smb (Q, E (Q, P) 
Q,rel (p_lq>P-Lq) 



^kl 



1 'T-gas 

* ijkl 

£ij=£, 



2 Q^/m^ 
dpf^p {p±q) 



Q 

2 



-Ski)e 

^kl 

3V^ 



-iQX/h 



Q 



Jkl 



rcl {p±q,P±q) - -77 



x/y (^rel(p^Q,P_LQ) - 
xSmb {Q,E{Q, ?)+£., 



Q 

2 
Q 

2 



kl 



E. 



(p_lq,P-Lq) 
Q,rel (pj_q,P_lq) 



Q 

2 
Q 

2 



Eh 



Q'^ /m, 
F ■ 

'-"1.3 

Q"^ jra. 



-Q 



Q 



'EliE,j,gj , 



(21) 



where X and P denote position and momentum operator of the massive test particle, and the scattering amplitudes fij 
appearing operator-valued describe inelastic scattering with a momentum transfer Q, between two channels differing 
in energy by £ij . One immediately checks that the diagonal matrix elements in the momentum representation of Eq. 
(I^Tj) do coincide with Eq. (PO)) . and furthermore that neglecting the internal degrees of freedom one comes back to 
the quantum linear Boltzmann equation [2^ . which together with the correct behavior under translations, granted 
by the very operator structure of the equation, provides a further argument for the assessment of the off-diagonal 
matrix elements. The step leading from Eq. (|20p to Eq. (PT|) . which corresponds to promote the classical momentum 
to the corresponding operator, similarly to what happens in standard quantization procedures, relies on the specific 
structure of Eq. (PO)) . and can also be applied when neglecting the internal degrees of freedom, in which case it 
leads to the correct version of the quantum linear Boltzmann equation, as confirmed by the independent derivations 
[igl . [2ll . [2^ . [30| . Of course the ultimate justification for Eq. (|21D relies on a microscopic derivation, which can be 
obtained similarly but with much lengthier calculations than in 22]. Due to the quite complicated expression it is 
worth introducing a more compact notation by defining the Lindblad operators 



-Q.P,£ = e*«-^/'"'L(p,P;Q,£:), 



(22) 



where 



L (p, P; Q,S) = ( "^el (P^Q' 



£ij=£ 



^ /m, 



-Q 



(p±q) J Smb (Q, E (Q, P) + £.,,)E,. 



thus writing Eq. (|2ip in the compact and manifestly Lindblad form 



dt 



j '^^ j ^Q,P,£e^Q,p,£ - \ {^Q.p,£^Q.P,£^ q) 



(23) 



(24) 



We will refer to Eq. (|24p or equivalently Eq. (PT|) as quantum Bloch-Boltzmann equation. 
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IV. REDUCED NON-MARKOVIAN DYNAMICS 



We now want to point out two different situations of physical relevance in which relying on Eq. (|24p one can 
obtain a description of non-Markovian behaviors, typically showing up in non exponential decay e.g. of coherences of 
the system under study. Despite focussing on a concrete class of physical systems, our analysis generally applies to 
the case of a bipartite quantum system interacting with a reservoir, provided all degrees of freedom of the bipartite 
system are involved in the interaction mechanism between system and reservoir, thus generating the entanglement 
which accounts for the memory effects. This provides a realization and clarification of the scheme envisaged in 
calling for extra fictitious unobserved degrees of freedom in order to lead to a Lindblad rate equation realizing in the 
Born approximation a generalized Lindblad structure. In particular our result goes beyond the Born approximation 
and displays the full-fledged generalized Lindblad structure Eq. ([5|), allowing for truly coupled equations for the 
different subcollections pa and considering both the case of a discrete and a continuous label a. 

As we discussed in Sect. U the non-Markovian behavior described via Eq. ^ arises when one goes over from 
Eq. p4)) to a semiclassical description, and a classical label characterizing the initial state cannot be resolved or 
accounted for in the final measurement. This semiclassical picture of the dynamics holds if the initial state of the 
system is prepared so that one of the two degrees of freedom is in a classical state, or if decoherence affects the two 
kind of degrees of freedom of the bipartite state on different time scales, so that e.g. the motional dynamics can be 
treated classically while the internal degrees of freedom still require a full quantum treatment. In this framework 
knowledge about the way in which the system is prepared usually naturally provides information about both parts of 
the bipartite system, while the final detection scheme is not necessarily fine enough to fully characterize the outgoing 
state. In different contexts it might also possibly arise as a consequence of superselection rules. 



A. Description of centre of mass degrees of freedom 



Let us consider first a situation in which we put our test particle, or equivalently a sufficiently dilute collection 
of such test particles so that they can be considered as non interacting, in a dense inert gas. The test particle will 
undergo many collisions quickly leading to a classical characterization of the motion of its centre of mass, so that 
only the diagonal matrix elements in the momentum representation of Eq. (j24p are left on a time scale set by the 
collisional decoherence mechanism, which leads us back to Eq. ([^0]) . which is also called Bloch-Boltzmann equation. 
Of course due to the complexity of the master equation Eq. (|24p such a behavior, though naturally expected on 
physical grounds and usually invoked in the literature on decoherence [sg'] , cannot be easily demonstrated in realistic 
situations. It has however been confirmed by means of Monte Carlo simulations, which also allow for estimates of 
the decoherence rates [STj . It is convenient to write the equation in a more compact way introducing the following 
(D-number rate operators 

M^l (P + Q;Q) = Se^^^s,, ^ / dp ^lp (p^q) Smb (Q, E (Q, P) + (25) 
X/.. frel(p^Q,P^Q) - ^ -f ^^Q,rel(p^Q,P^Q) + | 









Q2/m 




Ski 







xfu [rel P.q) - f + rel P.q) + f + 

which provide the rates for scattering from P to P + Q, including the dependence on the indexes for internal degrees 
of freedom and the transition energy £ij . In the limit of an infinitely massive test particle these rate coefficients can 
be checked to go over to those derived in [26*1 for the case of an immobile system. Exploiting the expression Eq. pS)) 
for the rate operators we can write Eq. ([20ll more compactly as 

^p(P)=^y'dQ M^UP;Q)E,,p{P-Q)El-^Mll{P + Q;Q)[ElE,,,piP)} . (26) 

ijkl 

Despite the enormous complexity of this integro-differential operator equation one can consider some simpli- 
fied situations, allowing to put into evidence non-Markovian behaviors arising for the reduced statistical operator 
p ~ j dPp{P) , only describing the internal degrees of freedom, once the dynamics of the various subcollections 
{/9(P)}pgj^3 is given by Eq. ([26]) . Let us consider the most simple conceivable situation, taking an internal (D^ space 
and only allowing for elastic scattering. We thus have M^j^ oc SijSki, and further restricting to forward scattering we 
can write 

Mg{P + Q-Q) = 6,,5ki5HQ)i^k{P). 
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parametrizing the rate operators by means of the functions S^ik (P). According to Eq. ((T2)) we denote by En = 

the maps between the same energy eigenstates, corresponding to the projectors on the two one-dimensional subspaces 

of C^ so that Eq. ^ now simphfies to 



d_ 
dt 



p{P) = (P) Eup{P)Ekk - ^ I E (^) p(-p) } • (27) 

ik \ i ) 



The dynamics of the single subcoUections p{P), can now be easily studied. Setting 

P,,(P) = {MP)\j) 

for the matrix elements of the collection {p (-P)}pg][{3 of matrices in (D^, one immediately sees that there is no 
dynamics for the populations, in that Pii{P) = 0, for j = 1,2, so that in particular integrating over the possible 
momentum dependence also pu = 0, for i — 1,2. The coherences of the subcoUections pi2{P) — P2i{P) are instead 
described by the equation 



-\s.ll{P)-\^22 {P)+il2 (P) 



P12{P), (28) 



which introducing what we might call, in the absence of better names, a momentum dependent friction coefficient 

S(P) = (P) + ^^2 (P) - ei2 (P) , (29) 

which in view of Eq. (I25p has a positive real part proportional to the averaged modulus of a difference of forward 
scattering amplitudes [26|], is easily solved by 

pi2(P,t) = e--(^)Vi2(P,0). 

Let us now consider an initial state of the form /9i2(P,0) = pi2{Q)pfj (P), corresponding to a preparation in which 
the test particle, or equivalently the dilute ensemble of non interacting test particles, is in a classical thermal state as 
far as centre of mass is concerned, and has a non vanishing initial value for the coherences of the internal degrees of 
freedom. The behavior in time of the off-diagonal matrix elements observed for the internal degrees of freedom only, 
not resolving the momentum of the considered test particle, is then given by 

pu{t) - A(t)pi2(0), 

with 

A(t) = J dPe--(^'V^(J')- (30) 

It is now immediately evident that behaviors utterly different from the usual Markovian exponential decay in time 
appear depending on the actual expression of the momentum dependent S (P) , the Markovian case obviously corre- 
sponding to a constant friction coefficient S (P) = rj. 

For the simple case S (P) = aP^ one immediately obtains a power law decay of the form 

= 7 (31) 

where we have set r = 1/aP^, indicating a natural reference time. Another simple expression of the friction coefficient 
leads instead to a stretched exponential. Considering in fact S (P) = 6/P^ one has to evaluate 

A(^) = ^/dPe-(-)*e--VP|, (32) 

so that exploiting the result jsl] 
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we obtain 



A(i) = l + (</r) 



.1/2 



-(t/r) 



1/2 



(33) 



describing a stretched exponential decay in time with a square root correction, where the reference time is now set by 
T — P^/ (46). These two simple choices for the friction coefficient S (P), amenable to an analytical treatment, have 
clearly shown the appearance of strongly non-Markovian behaviors for the operator p (t) . The considered example 
is obviously quite simplified and does not describe in a realistic way all possible aspects of the dynamics, e.g. the 
redistribution of population in the internal degrees of freedom. It allows however to easily grasp some non-Markovian 
aspects of the generalized Lindblad structure given by Eq. ([5]), of which Eq. (I26p provides a simple example, even 
though with a continuous index. 



B. Description of internal degrees of freedom 



We now focus on a quite different situation, in which we study the dynamics of our test particle when flying through 
an interferometer for massive particles, e.g. of the Mach-Zender type, as recently realized also for the quantitative 
study of decoherence [39j . In such a case the initial preparation is engineered so as to ensure a coherent superposition 
of states of the motional degrees of freedom of the system while, in the absence of a further selection in the prepared 
state, the internal degrees of freedom can be described by a classical distribution, corresponding to a partially diagonal 
statistical operator. The diagonal matrix elements of Eq. ([^^ with respect to the internal degrees of freedom naturally 
lead to coupled master equations for the collection {pr}^. of trace class operators in {^'^) : defined according to 
Pr ~ (''Ipk)- III order to keep a compact notation we introduce the rate operators 



R''Mp,Q) = /r, rel(p^Q,P^Q 



Q 



iQX/h 



^Q,rel(p,^,P,c,) + ^ 
Smb {Q,E{Q,P)+£r,), 



-Q 



(34) 



operator-valued on > so that the master equations for the subcoUections {pr}^ explicitly exhibit the generalized 

Lindblad structure given by Eq. ([5|) , with the additional appearance of integrals over continuous indexes 



dt 



Pr = 



J2 [dQ f dp \r^^ (p, Q) p,R^^ (p, Q)^ - i {r-'"'- (p, Q)^ R^^ (p, Q) , p,,} 



(35) 



The set of master equations given by Eq. (|35p do provide a non-Markovian dynamics for the statistical operator 
p — X]r=i observed at the outcome of the experiment, e.g. to determine the visibility of the interference fringes, 
when the detection scheme cannot resolve the state of the internal degrees of freedom. 

In a typical experimental situation for the study of coUisional decoherence one can safely neglect the dependence 
on the momentum operator P in the rate operators defined by Eq. I|34p. replacing it by the classical value of the 
momentum of the incoming test particle, due the fact that on the decoherence time scale the dissipative dynamics 
of the momentum does not play a role (2^ . l40l . l4lj | . This brings in an important simplification in Eq. (|35p , which 
can now be written by putting into evidence the unitary operators e^Q'^^^, describing the momentum kicks causing 
decoherence of the centre of mass, and (D-number positive collision rates Xjr (Q), which depend on the internal state 
of the test particle, thus obtaining 



dt 



Pr 



^ J2 dQ IXrj (Q) e^'^ ^'/^pje-'^ ''/^ - A,, (Q) pr 



(36) 



For the case in which the collisions do not lead to transitions between different internal states, so that Aj> (Q) = 
SjrXr (Q)i one comes to the following master equation describing a dynamics determined by momentum kicks of 
amount Q, taking place with a probability density Vr [Q) which depends on the internal state of the test particle 



dt 



Pr 



K / dQVr iQ) 



'■x/h 



Pre 



In Eq. ([37| the probability density Vr (Q) is defined according to 

A. (Q) 



VriQ) 



Ar 



'■X/h 



Pr 



(37) 



(38) 
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with Ar = J dQ\r (Q) the total scattering rate for particles with internal state r. The master equations Eq. (|37|) are 
easily solved in the position representation according to 

{x\pr{t)\y) = c-^'^[i-^-(--^)l*(a^|p,(0)|y), (39) 

with 

the Fourier transform of the probability density given by Eq. ([38|) . that is to say its characteristic function [i^ . 
In particular the prefactor is the characteristic function of a compound Poisson process, composed according to the 
probability density Vr (Q)- This equation describes a quite general physical situation in which one has a sequence 
of interaction events between system and environment, distributed in time according to a Poisson distribution, each 
one characterized by a random momentum transfer, drawn according to a certain probability density fixed by the 
microphysical interaction mechanism (43l . . At variance with a simple Poisson process the momentum transfer is 
not deterministically fixed to be the same in each collision, but is a random variable depending on the details of the 
collision. 

We now look at the dynamics of the matrix elements of the whole statistical operator p = Pr i responsible for the 
description of the measurement outcomes at the output of the interferometer. To do this we consider an initial state 
of the form pr (0) — Prp{0), with Pr 1^ and ^j,Pr = 1 , corresponding to a preparation in which the populations 
in the internal states are distributed according to certain weights Pr, while the statistical operator p{0) in {^^) 
characterizes the quantum state of the centre of mass. We thus come to the following expression for the solution of 
the non-Markovian set of master equations given by Eq. ([57]) 



{x\p{t)\y) = ^p,e-A'-[i-*-(-2')l*(a;|p(0)|y), (40) 

r 

where the multiplicative prefactor determines the weight of the matrix elements in the position representation, both 
with elapsing time, and as a function of the distance \x — y\. Such an expression provides information on the loss of 
coherence responsible for reduction in the visibility of the interference fringes. If only one of the weights Pr is different 
from zero, and therefore equal to one, one falls back to the usual Markovian exponential decay in time, possibly with 
a modulation in the spatial dependence. In particular the coherence of the quantum state over spatially separated 
points depends on the details of the functions (pr {x — y), given by the Fourier transform of the probability density of 
momentum kicks in the scattering events. In a typical situation such functions quickly go to zero in their dependence 
on the distance |a; — y|, so that e.g. in the Markovian case one is simply left with a constant exponential loss of 
visibility [3^. With the present more general initial state even for (pr (x — y) ~ one has a non trivial structure 
describing a decay of coherence other than exponential 

{x\p{t)\y) 0, ^it){x\p{0)\y}, (41) 

where the function 



^(t) = ^p.e-^-*, (42) 

r 

is the survival probability of a multiexponential distribution, i.e. the probability to have no event up to a time t 
for such a distribution '42]. Depending on the weights Pr and the rates A^, not only simple deviations from the 
exponential law can appear, but also utterly different behaviors. To clarify this point let us consider for the sake of 
example a geometric distribution of weights, with ratio po = e^", a G IR+, so that 

Pr = (1 -po)Po, 

and a geometric progression of rates 

K = Ao75, 

with ratio 70 — e^^, b G R+, and the reference rate Aq as scale factor. The survival probability then reads 



r 
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which due to the relation [i^ 

«'(Ao<) = l[vl/(i)-(i-po)e-Ao*] 
Po 

exhibits at long times a power law decay 

f (Aoi) 'i^ (43) 

This most simple example allowing for an analytical treatment already shows the rich variety of non-Markovian 
behaviors which might arise when one uses the generalized Lindblad structure given by Eq. (|35p . which provides a 
further example of the general result Eq. ([5]) for the case of a sum over a discrete index. In particular Eq. pS]) 
describes how the dependence of the scattering events on the internal structure of the test particle affects the loss 
of coherence in position space, which in turn determines the reduction of visibility in an interferometric experiment. 
As it appears from Eq. (j43p this can lead to very strong deviations from the exponential decay, such as power law 
behaviors. 



V. CONCLUSIONS 



In this paper we have considered a class of non-Markovian behaviors, arising when dealing with a bipartite quantum 
system interacting with a reservoir. The concrete bipartite system considered was given by a massive test particle, for 
which both internal and centre of mass degrees of freedom have been taken into account. The reservoir was assumed 
as a structureless gas, affecting our test particle through collisions whose microscopic characterization depends on 
both its motional and internal state. As a starting point we have derived in Sect. Illll a new quantum master equation 
describing such dynamics in a non perturbative way, expressed by Eq. ([M)) . which can also be termed quantum 
Bloch-Boltzmann equation in that it describes at the quantum level both kinds of degrees of freedom. When due 
to decohercnce or features of the initial preparation one of the two degrees of freedom is to be described classically, 
one obtains from the quantum Bloch-Boltzmann equation two examples of a generalized Lindblad structure recently 
considered for the description of non-Markovian dynamics 0, Q • Such a generalized Lindblad expression has been 
outlined in Sect. |TT1 clarifying its mathematical structure and physical motivation. For the case at hand non-Markovian 
effects, leading to decay behaviors of coherences of the system given by stretched exponentials or power laws instead 
of simple exponentials, appear when the degrees of freedom allowing for a classical description are not resolved in 
the final measurement, only focussing on the quantum degrees of freedom. This provides a concrete realization of a 
proposed mechanism for the appearance of such generalized Lindblad structures [isi . [l6j , further clarifying the origin 
of the non-Markovian behaviors. These behaviors have been spelled out in Sect. IIVI focussing on the dynamics of the 
internal state of a test particle interacting with an inert gas, as well as on loss of coherence of a massive particle flying 
through an interferometer where it interacts with a background gas. It is to be stressed that in the physical examples 
considered in Sec. IIV Al and Sec. IIVBI one has to deal with decoupled subcollections of statistical operators, so that 
the non-Markovian features arise from the average over the classical index in the initial condition. The generalized 
Lindblad structure given by Eq. (O also allows to consider coupled equations for the different subcollections, and for 
these situations one naturally expects a much more complicated non-Markovian dynamics. 

It immediately appears that the outlined scheme leading to a class of non-Markovian evolutions generally applies in 
the presence of the interaction of a bipartite quantum system with a quantum environment, when one of the quantum 
labels of the system becomes classical and can be averaged over. More generally such a class of non-Markovian 
evolutions appear in the presence of a classical degree of freedom, described by means of some discrete or continuous 
label, which is involved in the characterization of the interaction between two quantum degrees of freedom, and 
is averaged over in order to give the relevant dynamics. This classical label might as well appear on the side of 
the environment, corresponding to so called structured reservoirs, or on the side of the system, as in the case of a 
bipartite system. The present work can naturally be extended to include an internal structure in the gas particles, 
which could also influence the scattering amplitude, introducing new channels. In particular a detailed analysis of the 
rate operators based on microphysical informations could pave the way to new interferometric experiments for the 
quantitative study of decoherence, exhibiting more general behaviors than exponential decay of visibility with elapsed 
interaction time. 
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